New family of models for incompressible quantum liquids in d > 2 
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Through Haldane's construction, the fractional quantum Hall states on a two-sphere was shown 
to be the ground states of one-dimensional SU(2) spin Hamiltonians. In this Letter we generalize 
this construction to obtain a new class of SU(N) spin Hamiltonians. These Hamiltonians describes 
center-of-mass-position conserving pair hopping fermions in space dimension d > 2. 



Gapped quantum many-body systems are stable states 
of matter because they are robust against weak pertur- 
bations. For fermions, the most common gapped system 
is the "band insulator" where an energy gap in the dis- 
persion relation separates filled and empty single-particle 
states. Here the insulating behavior is caused by the 
Pauli exclusion principle. We refer to a system as a 
"many-body insulator" if its energy gap is caused by 
many-body interaction rather than one-particle disper- 
sion relation. It is widely believed that such a many-body 
gap can exist when the occupation number (i.e., the aver- 
aged number of particle per spin per unit cell) is a frac- 
tion. The Mott insulator, where the energy gap is due to 
the inter-particle repulsion, is an example of many-body 
insulator. 

Contrary to the common belief, it is very difficult to 
find true many-body insulators. In most cases an energy 
gap at fractional occupation number is accompanied by 
the breaking of translation symmetry. After symmetry 
breaking, the unit cell is enlarged such that the new occu- 
pation number is an integer. In a recent work Oshikawa 
argued that the existence of energy gap at fractional oc- 
cupation requires the ground state to be degenerate [1]. It 
happens that under usual circumstances such degeneracy 
is caused by translation symmetry breaking. 

The fractional quantum Hall state is a quantum liq- 
uid (hence no breaking of translation symmetry) with 
an energy gap. On the surface it is not clear what does 
it have to do with the many-body insulator discussed 
above. In Ref.[2-4] it is shown that, when placed on a 
torus, both abelian and non-abelian quantum Hall liquids 
can be mapped to many-body insulators on a one di- 
mensional ring of lattice sites. When both dimensions of 
the torus are much bigger than the magnetic length, the 
energy gap is generated by a long range center-of-mass- 
position conserving hopping, rather than density-density 
interaction. 

Partly motivated by Anderson's proposal of spin 
liquid [5], the question of whether a many-body insula- 
tor can exist without symmetry breaking in spatial di- 
mension d > 2 has attracted a lot of interests. In this 
Letter we give this question an affirmative answer by ex- 
plicitly constructing a new class of solvable lattice mod- 
els that exhibit incompressible quantum liquid ground 
states. Since our construction is a generalization of Hal- 



dane's work on the pseudopotential Hamiltonian for the 
fractional quantum Hall effect [6], we shall begin by briefly 
review it. 

If a magnetic monopole of strength 2S (5 is a multi- 
ple of 1/2) is placed at the center of a two-sphere, the 
kinetic energy spectrum of a particle confined to move 
on the sphere is given by Ek oc (S + k)(S + k + I), 
where k is an non-negative integer. The kth "Landau 
level" is 2(5 + k) + 1-fold degenerate. The degeneracy 
of the lowest Landau level, k — 0, is exactly the dimen- 
sion of a spin-S SU(2) multiplet. Thus the Hilbert space 
of N spin polarized electrons in the lowest Landau level 
is the same as the exchange-antisymmetric sub-Hilbert 
space of N such SU(2) spins. Haldane's pseudopotential 
Hamiltonian (which has the spherical version of Laugh- 
lin's v = 1/ra wavefunction as the ground state) is given 
by 
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Here i,j — 1, .., N are the spin labels, -P? 5-9 projects the 
product states of spin i and j onto the total spin 25 — q 
multiplet, and n q > are parameters (as a result Eq. (1) 
is positive-definite). In Eq.(l) the q-sum is restricted to 
odd integers because the restriction of the Hilbert space 
to the total antisymmetric subspace. It can be shown 
that when N — 1 = 2S/m = p the Hamiltonian in Eq.(l) 
has an unique ground state described by the following 
spin coherent-state wavcfunction(thc spherical version of 
Laughlin's wavefunction) 
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Upon expanding the determinant, Eq.(2) can be written 
as a linear combination of JI^Li V (rnp)\ / njlkjlu" 3 v^ 3 
where nj+kj = mp. Each term in this linear combination 
is a direct product of N, spin-S*, states. An important 
property of the wavefunction in Eq.(2) is that the 
highest total spin for any pair is 25 — m. Consequently 
Eq.(2) is a zero energy state of Eq. (1). Moreover, if we 
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view the 2S + 1 different S z states as the 2S 
orbitals of a one-dimensional lattice 



f local 
and write Pff 9 = 

A= — {2S—q) Z-^rri2 — — S Z-^m-i — — S u 5,mi .5,i-mi u S,m2,5,(-ni2 

m 2 >< mi, I — mi\, Eq.(l) becomes a center-of-mass 
conserving pair hopping Hamiltonian.[2, 3] In the above 
^Sm'si-m, ^ S tne SU(2) Clebsch-Gordon coefficient. 
The role of center-of-mass position conservation in 
producing true many-body insulators was discussed in 
Ref.[2, 3]. 

In the following we generalize Haldane's construction 
to SU(3) spins. The reason for doing so is SU(3), a rank 
two Lie group, has multiplcts isomorphic to two dimen- 
sional lattices. Hence it allows us the possibility of con- 
structing Hamiltonians for many-body insulator in d = 2. 
The irreducible representations of SU(3) are labelled by 
two integers (p,q). In the following we shall focus on 
the the multiplets (k,0). The reason is because these 
multiplets are the only ones whose weight space is an 
array of non- degenerate, i.e., non-duplicated, points (see 
Fig.(l)). In the following we consider SU(3) spins each in 
the (mp, 0) representation. The dimension of the (mp, 0) 
representation is d(mp) = (mp+l)(mp+2)/2. For reason 
that shall become clear later, we shall choose the num- 
ber of spins so that N = d(p) = (p + l)(p + 2)/2. Under 
such condition the filling factor, / = d(p)/d(mp), is 1/m 2 
in the thermodynamic (p — > oo) limit. As earlier, we 
will constrain the N-spin Hilbcrt space to be exchange- 
antisymmetric to mimic the fermion statistics. 

The spin Hamiltonian we construct is a generalization 
of Eq. (1), and is given by 
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Here the operator p@ mp ~ 2q ' q ) operates on the direct 
product states of two spins i and j, and projects them 
onto the (2mp — 2q,q) multiplet, and n q > 0. For sim- 
plicity in the rest of the Letter we shall set m = 3 and 
p = odd integer. In this case Eq.(3) becomes 
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FIG. 1: The weight space of (6,0). 



by the following SU(3) coherent-state wavefunction 
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First we prove Eq.(5) is a ground state. Let us focus on 
the dependence of Eq. (5) on the variables of any chosen 
pair of spin i and j. For each of the Slater determinant in 
Eq. (5) the highest total SU(3) weight of these two spins 
is (2p — 2, 1) because of antisymmetry As a result when 
we multiply three determinant together the highest total 
SU(3) weight for spin i and j is (6p — 6, 3). Consequently 
Eq.(5) is an zero-energy eigenstate of the positive-definite 
Hamiltonian in Eq.(4). Thus we have found a ground 
state. 

As to the uniqueness let us start with the simplest case 
of p = 1 where the single spin Hilbert space is d(3) = 10 
dimensional and there are d(l) = 3 spins. The most 
general many-body wave function is given by 
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where ^i' 2,3 = u 



, Vj , Wj . The requirement that this 
wavefunction lies in the direct product of (3, 0) <S> (3, 0) (g> 
(3, 0) demands the coefficient C to be invariant when the 
three indices of any chosen particle are permuted. In ad- 
dition, the antisymmetry constraint restricts C to change 
sign upon the exchange of particle labels. Next, let us 
pick any pair of spins i and j and examine the dependence 
of Eq. (6) on their variables. The possible total SU(3) 
weights of these two spins that are consistent with anti- 
symmetry are given by (3, 0) (8) (3, 0) = (4, 1) © (0, 3). In 
order for the wavefunction to be annihilated by Eq. (4), 
it must lie entirely in (0,3). Since there are only three 
particles it is not hard to show that there is a unique C 
so that the above condition holds for all pair (i, j): 

C(ai, 0.2, C*3; Pi, 02, 03', 7l7 72, 73) OC eQ 1 /3i7iea 2 /3 2 7 2 ea 3 /3 3 73- 

Substitute the above equation into Eq. (6) gives \ ~ *3- 
Now, we consider p = 3, where the single-particle 
Hilbert space is 55 dimensional and and there are c?(3) = 
10 particles. Analogous to Eq. (6) the most general 10- 
particle wavefunction is given by 
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(7) 



{a 3 „ = l} 



j=l n=l 



We shall spend much of the rest of the Letter to prove 
that Eq.(4) has a unique singlet ground state described 



Here j is the particle label, n = 1, ...,9 labels the nine 
fundamental SU(3) spinors (1,0) that make up (9,0). 
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The symmetry properties of C are the same as be- 
fore. In this case the total SU(3) weight of any two 
spins that are consistent with antisymmetry are given 
(9, 0) ® (9, 0) = (16, 1) © (12, 3) © (8, 5) © (4, 7) © (0, 9). 
The condition of being annihilated by the Hamiltonian 
requires the ground state wavefunction to lie entirely in 
(12, 3) ©(8, 5) ©(4, 7) ©(0,9), i.e., 

{P^+P$*+P^+P$»)X = X, V(U). (8) 

Eq. (8) implies that among the 9 indices spin i and j 
each possesses, there must be at least 3 pairs (a pair con- 
tains one index from each particle) such that C — > — C 
upon exchanging the indices within each pair. In addi- 
tion to be consistent with exchange antisymmetry, the 
number of such pairs also must be odd. Since C is in- 
variant when the nine indices of any particle are per- 
muted, we can perform permutations so that in each 
triplet {a il7 ai2,a i3 ), (an, a i5 , a i6 ),(a i7 , a i8 , a i9 ) of par- 
ticle i and (ctji, ctj2, aj 3 ), (aji, ctjs, ajo),(aj 7 , ctjg, ajg) of 
particle j there is an odd number of antisymmetric in- 
dices. Under such circumstance C changes sign upon 
independent exchanges of the triplets, i.e., 

C — ► -C upon (aii,aj2,a»3) ^ ( a ji, a j2 , a j3 ) 
(an, a»5, a i6 ) <-► (0,4, a j5 , a j6 ) 
(a i7 ,a i8 ,a i9 ) (a j7 ,a js ,a j9 ). (9) 

In Eq. (9) (...) (...) denotes the exchange of whole 
group of indices. If Eq. (9) can be made true simultane- 
ously for all pairs i and j then ^3 is the unique solution 
of Eq. (8). This is proven as follows. 

Let us focus on the dependence of C on the first in- 
dex triplet of all particle. For each triplet of indices, 
say {a a ,ai2,a i3 ), there are (3 + 2)!/(3!2!) = d(3) = 10 
incquivalcnt combinations. We can interpret each com- 
bination as a single-particle quantum state and C as the 
wavefunction for 10 particles to occupy these states. The 
first line of Eq. (9) allows us to interpret C as the wave- 
function for fermions. For N = d(3) = 10, i.e., when 
the fermion number is the same as the number of single- 
particle state, there is an unique wavefunction satisfying 
the antisymmetric requirement, namely, 

C({an,ai2,ais}...) ~ e{(anai2a i3 )} (10) 

where e{...} is the rank 10 total antisymmetric tensor 
with respect to the exchange of index-triplets. Sim- 
ilar argument can be made to {(an, a^, a^)} an d 
{(a i7 , a i8 , a i9 )}, and lead to 

C({an, ...,a i9 }) ~ e{(aiiai2a a )}e{(a iA a i?> ai & )} 
x e{(a l7 a l8 a l9 )}. (11) 

Substitute Eq. (11) into Eq. (7) we obtain \ ~ *3- 

Now, we shall prove that Eq.(9) can indeed be made 
true for all pairs i and j simultaneously Let us assume 



that there exists a ground state solution whose C does 
not satisfy Eq.(9) for pair (k,l). This means there must 
be at least one triplet exchange, let say {aki, ak2, a k3 } 
{an, ai2, ai 3 }, for which C does not transform according 
to Eq. (9). However, since the wavefunction still has 
to satisfy Eq. (8) for (k,l), we should be able to write 
C = C3 + C5 + C 7 + Cg , where C q is the component of C 
that is odd with respect to exchange of exactly q pair of 
indices between particle k and I and even with respect to 
the exchange of the rest. Now let us consider the effect 
of {a kl ,a k2 ,a k3 } <-> {an, an, ai 3 } on C. Under such 
operation C q can either change sign or stay invariant de- 
pending on whether an odd or even number (out of q) an- 
tisymmetric indices are contained in the specified triplets. 
In other words upon {a k i, a k2 , 0^3} {an, an, ao} we 
have 



C -> n 3 C 3 + r) 5 C 5 + r] 7 C 7 + r] 9 C 9 , 



(12) 



where r\ q = ±1. Since Eq.(9) is not satisfied, 773,5,7,9 must 
not simultaneously be —1. Now consider a new C 



C' = 



1 r 



C - r] 3 C 3 - rj 5 C 5 - r] 7 C 7 - t] 9 C 9 



(13) 



It is obvious that upon {a k i, a k2 , a k3 } <-> {an, an, ai 3 } 
C — > —C. Moreover by construction C only contains 
those C q whose rj q = —1. Now use C as the starting C 
and repeat the above operation until we reach a final C 
for which Eq.(9) holds for all triplet exchanges and for 
all (i, j). Since at each stage of obtaining C certain C q 
are projected out, there must be missing q components in 
the final C. However we have already proven that any C 
that satisfy Eq.(9) for all (i,j) pair must satisfy Eq. (11). 
However Eq. (11) contains all four q components for all 
pair (i, j). Consequently we have reached a contradiction. 
Therefore it must be possible to make Eq.(9) hold true 
for all pairs (i,j) for any ground state solution satisfying 
Eq. (8). 

Although we have chosen p = 3 and m = 3 in the 
above discussion, it should be clear that our proof can be 
generalized to any odd p , any to. Thus we have proven 
that Eq.(5) is the unique ground state of Eq.(4). 

It is straightforward to prove that Eq.(3) is a center- 
of-mass position conserving pair hopping model, i.e., 



j,L,L 3 l,l 3 k,k. 



<3,L,L 3 rpj,L,L 3 
Ma r k,k 3 



The central steps are 1) viewing the weight space of 
(mp, 0) as a triangular lattice, and 2) decomposing the 
two spin states in Eq. (3) as linear combination of prod- 
ucts of single spin state. Due to space limitation, the 
result will be published elsewhere. Using the explicit ex- 
pression for the F's in Eq. (14) (lengthy hence is omitted 
here) we can estimate the hopping range. In Fig. (2) we 
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plot the hopping range versus p for p — 30 — > p = 600. 
While the linear dimension of the lattice scales as p, the 
hopping range scales as hence the hopping is long- 
ranged as the SU(2) case. 
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FIG. 2: A log- log plot of the hopping range and p from 30 to 
600. The straight line is the best fit. The vertical axis is the 
log of the hopping range and the horizontal one is log p. The 
hopping range scales as 

Finally, we demonstrate the presence of an excitation 
gap within the single mode approximation (SMA)[7]. 
Analogous to the SU(2) case it is possible to view the 
(k, 0) SU (3) multiplet as the "lowest Landau level" (LLL) 
of a particle running on CP 2 under the action of a U(l) 
background magnetic field [8]. If we parameterize the fun- 
damental SU(3) spinor as (l,zi,z 2 ) /y/X + \zi\ 2 + |z 2 | 2 , 
where Zi = Xi +iyi, and take the flat-space limit (i.e., re- 
strict |zi. 2 | << 1) the single-particle orbitals in the LLL 
become 
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where are non-negative integers. We recognize that 
the above result is the product of two LLL wavefunctions 
in two space dimensions. Thus in the flat-space limit, the 
LLL in CP 2 becomes the direct product of the LLLs in 
two quantum Hall planes. This reduction allows us to 
perform the SMA calculation pretty much in parallel to 
that for the ordinary quantum Hall effect. [7, 9] in the 
following we summarize the results. Within SMA the 
excitation energy is given by 



A(k) = /(k)/ s (k). 



(15) 



Here k = (k\, fc 2 ) where fci j2 are the complex wave vectors 
associated with the two quantum Hall planes, and /(k) 
and s(k) are given by 

f(k) = (l/N)<^ m \[pl[V, Pk }}\^ m > 

s(k) = (1/N) < * m |/^ / o k |* ro > . (16) 

In the above equation pk and V are the the density op- 
erator and the inter-particle potential projected onto the 



LLL. Straightforward calculation gives, 

1 X — / 1 i . . qk kq , 



/(k) = -5>(|q|)(e¥- 
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+s(k + q)(e ■ -eT)], 
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(17) 



where v(|q|) is the Fourier transformation of the poten- 
tial, which is required to be positive indicating the re- 
pulsive interaction to ensure the excitation energy to be 
positive. On the other hand, s(k) can be related to the 
radial distribution function g(f) by 



S (k) 



d 4 re 



—ik-r 



\g(f)-l]+p(2ir)W(k)(18) 



where p is the average density and k=(He(ki), Im(fci), 
Re(fc2), Im(fc 2 )). After some algebra, for small |k| 
it can be shown that A(k) = (a|fci| 4 + 6|fci| 2 |fc 2 | 2 + 
a|fc 2 | 4 )/(c|/ci| 4 + d|/fci| 2 |/fc 2 | 2 + c|fc 2 | 4 ), which remains fi- 
nite as k approaches to zero in any direction. 

In summary, we have constructed a two dimensional 
center-of-mass conserving pair hopping model which ex- 
hibit incompressible quantum liquid ground state. Al- 
though throughout the Letter we have focused on SU(3) 
whose weight space is two dimensional, our construction 
can easily be generalized to SU(N) giving rise to models 
for incompressible quantum liquid in higher dimensions. 
Particularly, the SU(4) model can be applied to the four- 
dimensional quantum Hall effect proposed by Zhang and 
Hu [10-12]. 
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